Polarization dependence of phonon influences 
in exciton-biexciton quantum dot systems 
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We report on a strong dependence of the phonon-induced damping of Rabi dynamics in an opti- 
cally driven exciton-biexciton quantum dot system on the polarization of the exciting pulse. While 
for a fixed pulse intensity the damping is maximal for linearly polarized excitation, it decreases 
with increasing ellipticity of the polarization. This finding is most remarkable considering that the 
carrier-phonon coupling is spin-independent. In addition to simulations based on a numerically 
exact real-time path integral approach, we present an analysis within a weak coupling theory that 
allows for analytical expressions for the pertinent damping rates. We demonstrate that an efficient 
coupling to the biexciton state is of central importance for the reported polarization dependen- 
cies. Further, we discuss influences of various system parameters and show that for finite biexciton 
binding energies Rabi scenarios differ qualitatively from the widely studied two-level dynamics. 



I. INTRODUCTION 



Investigations devoted to the carrier dynamics of op- 
tically driven quantum dots (QDs) with or without an 
additional coupling to a cavity constitute a highly ac- 
tive field of physics* 1 - - — In the limit of strong electronic 
confinement and for excitations with arbitrarily polar- 
ized laser pulses, the electronic structure of a QD can 
be modeled as a four-level system, consisting of the 
ground state (no electron- hole pair), two energetically 
almost degenerate single exciton states with orthogo- 
nal spin polarizations and the biexciton state. Such 
exciton-biexciton systems are attractive for many innova- 
tive technologies, including the realization of sources for 
on-demand entangled or squeezed photons^— or appli- 
cations in the field of quantum information theory* 1 ^— 
The major obstacle for coherent manipulations is deco- 
herence, and much effort has been devoted to character- 
ize or mimimizc decohcrence processes in optically driven 
QDs4"— i 11 ' 19 " — Recent experimental studies of excitonic 
Rabi rotations 4 ^ identified the pure dephasing coupling 
of carrier states to bulk like acoustic phonons as the prin- 
cipal source of dephasing in self-assembled QDs and con- 
firmed theoretically predicted results for phonon-couplcd 
two-level systems, such as a renormalization of the Rabi 
frequency 22 ' 24 Moreover, first evidence for the expected 
undamping of Rabi rotations at high pulse area o 3 ' 28 has 
been shown by observing a roll-off behavior in the oscil- 
lation amplitude^ 

Obviously, the dissipative dynamics of a four-level sys- 
tem accounting not only for single exciton but also for 
biexciton excitations can develop much more manifold 
than that of a two-level system and is in view of possible 
applications of special interest. Most of the existing lit- 
erature on the phonon impact on the carrier dynamics of 
optically driven QDs, however, deals with electronic two- 
level systems and comparatively little is known about the 
more complex exciton-biexciton dynamics. In Refs. 
and analytic expressions have been derived for the 
limit of ultrashort excitations. Two-photon Rabi oscil- 



lations between the ground and the biexciton state have 
been studied experimentally in Ref. [3l| and the phonon 
influence on this scenario has been analyzed theoreti- 
cally within a perturbative approach in Ref. [32l . Very 
recently, we have demonstrated that a numerically ex- 
act real-time path integral approach previously used for 
two-level systems ^ 28 ' 33 can be also implemented for four- 
level systems paving the way towards investigations in 
the regimes of high temperatures, strong carrier-phonon 
couplings and long times for arbitrary optical drivings^ 4 - 
In that work dark superpositions were identified as a gen- 
eral mechanism that strongly modifies the quantum dissi- 
pative relaxation, in particular at long times, and various 
related dynamical features have been discussed, including 
a crossover between two qualitatively different relaxation 
scenarios^ 4 - 

Here, we analyze in detail the impact of acoustic 
phonons on the dynamics on shorter time-scales than 
in Ref. concentrating on Rabi scenarios and signals 
most convenient to be studied in experiments. In partic- 
ular, we study the dynamics for different intensities and 
polarizations of the driving laser field. As the exciton- 
phonon coupling is independent of the spin, it might be 
suspected, that the phonon influence does not depend on 
the polarization of the optical excitation which creates 
well defined superpositions of carrier spins. Nevertheless, 
we shall demonstrate that in contrast to this expecta- 
tion, the phonon-induced damping docs strongly depend 
on the polarization of the exciting pulse. Keeping the 
pulse intensity fixed, the damping is most pronounced 
for linear polarization, while for elliptical polarizations it 
becomes the more reduced the closer the polarization is 
to the circular limit. 

The paper is organized as follows. In Sec. |IT]we intro- 
duce the model, we present a short summary of the path 
integral approach used and comment on the challenges 
that are faced when accounting for four electronic levels 
within this numerically exact formalism. In addition, we 
introduce a weak coupling theory that allows for approx- 
imate but explicit expressions for the pertinent damping 
rates. This perturbative approach is described in more 
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detail in the appendix and most clearly reveals the ori- 
gin of the polarization dependence of the phonon-induced 
damping. The latter is presented in Sec. IIII[ where we 
also study in detail the dynamics for linearly polarized 
excitations and discuss the influence of pulse intensity, 
temperature and magnitude of the biexciton binding en- 
ergy. Special emphasis is put on comparisons with results 
for two- level systems. Finally, we summarize our results 
in Sec. HVl 



mechanism^!. With 6 q denoting the creation operator 
of a LA phonon with wave vector q and energy fiw q , -ff p h 
and -f/dot-ph are then given by 



P h 



Wq&q&q , 



#dot-ph = ( 7 q & 1 + ^ 6 q) n "\ V ){ V \ 



(4) 



(5) 



II. THEORY 
A. Model 

We consider a strongly confined GaAs QD with spin 
degenerate electronic particle states coupled to a contin- 
uum of acoustic phonons and driven by laser light. The 
corresponding Hamiltonian can be written as 

H = Hdot + H p h + -ffdot-ph + -f/dot-light , (1) 

where f/dot describes the electronic structure of the QD, 
.ffph is the free phonon Hamiltonian, -f/dot-ph represents 
the carrier-phonon and Hdot-iight the carrier-light cou- 
pling. As in the strong confinement limit the electronic 
single particle states are energetically well separated, we 
can safely neglect a Coulomb-induced mixing of states 
with different single particle energies and concentrate on 
the topmost valence and lowest conduction band states. 
Neglecting excitonic states that are not coupled to other 
electronic states via the carrier-light interaction, our sys- 
tem comprises four states. Besides the ground state |G) 
without electron-hole pairs we account for the two single 
exciton states |+) and |— ) with a total angular momen- 
tum of ±1, and the biexciton state |B). Then, i?dot reads 



The dispersion is assumed to be linear, cu q 



|q|, where 



ffdot = 7*n(|+)(+| 



-)(-|) + (2«1 

l-X+l) ■ 



A)|B)(B| (2) 



We take the unexcited ground state as the zero of energy, 
Htt defines the ground state exciton transition energy, A 
denotes the biexciton binding energy resulting from elec- 
trostatic Coulomb interactions and V cx the electron-hole 
exchange interaction. For simplicity, we consider spher- 
ical dots with envelope functions given by the ground 
state solution of a harmonic potential, i. e. 



■<Pe{h)(r) = 



1 



e(h) 



exp 



e(h) , 



(3) 



where a e fh\ denote the localization lengths of electrons 
and holes, respectively. 

To model the carrier-phonon interaction we concen- 
trate on pure dephasing coupling s 19 ' 35 of the electronic 
states to bulk like longitudinal acoustic (LA) phonons 
via the deformation potential that for strongly confined 
GaAs QDs has been identified as the dominant dephasing 



v c = 5110 m/s is the sound velocity, n v is the number of 
electron-hole pairs present in the state \v) (i. e., uq = 0, 
n+ = rt_ = 1 and ub = 2) and 7 q denote the exciton- 
phonon coupling constants. We stress, that the carrier- 
phonon coupling as given in Eq. ([5]) is independent of the 
carrier spin. Besides, it should be noted that the factor- 
ization in Eq. ([5]) into a procuct of 7 q and n v is valid for 
strongly confined QDs, where it is justified to treat exci- 
tons and biexcitons as product states representing uncor- 
relatcd electron-hole pairs. However, this factorization 
does not hold in general. For QDs far from the strong 
confinement limit, c. g., interface QDs, measurements 3 ^ 
indicate strong deviations from Eq. ([5]). Details on the 
phonon coupling can be found in Ref. Il9l . 

Finally, the coupling to the light field reads within the 
common dipole and rotating wave approximation 

ffdot, light = -h[F a+ (t)(\+)(G\ + |B)H) (6) 
+ J F CT _(t)(|-)(G| + |B)(+|) + h.c], 



where F a± = (t)M/h with E^ ; (t) being the positive 
frequency part of the o± circularly polarized component 
of the light field and M denoting the ground state ex- 
citon transition dipole moment, that we assume to be 
identical to the exciton biexciton transition dipole mo- 
ment. Taking the laser frequency to be in resonance 
with the polaron shifted exciton transition frequency 
H = fl - E q 7qMj, we set F a± (t) = f a± (i) exp(-»nt)/2 
with real envelope functions f a± (t) , that we will refer to 
as <j± circularly polarized field strengths. The total field 

strength is given by f(t) = , / (f ) + (t) and defines 



(+)/ 



(7) 



the pulse area a of an applied laser pulse via 



f(t)dt. 



For linearly polarized laser pulses, the total pulse inten- 
sity is equally distributed between both circularly polar- 
ized components, i. e., f a+ — f a _ —f/y/2, and both spin 
degenerate excitons experience the same dynamics. For 
circularly polarized pulses with / = f a± and f a = 0, 
the system reduces for V cyi = to an effective two-level 
system. All other polarizations we will refer to as ellipti- 
cal. Note, that the definition of the pulse area in Eq. ([7]) 
is such, that for circular polarizations (i.e., a two- level 
model), a pulse area of a = it inverts the system^ 
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B. Path integral approach 

To study the combined carricr-phonon dynamics, we 
use a real-time path integral approach, that for the four- 
level system has first been applied in Ref. O This ap- 
proach is numerically exact in the sense that there are no 
approximations within the model given above, and thus 
we face no errors beyond well controllable discretization 
errors. The reduced electronic density matrix is calcu- 
lated by taking the Liouville von Neumann equation for 
the full density matrix and tracing over the phonon de- 
grees of freedom, where the time evolution operator is 
represented as a path integral. Tracing out the phonon 
degrees of freedom introduces a memory and leads to a 
non-Markovian dynamics for the reduced density matrix, 
that can be fully accounted for. A complete description 
of our algorithm can be found in Ref. l38l where we give 
not only a detailed derivation, but highlight also specifics 
of our implementation that are related to the superohmic 
character of the carrier-phonon coupling and that are 
not encountered in implementations of path integrals for 
studies with ohmic or subohmic coupling types^r— 

The numerically complete treatment within the path 
integral approach rapidly becomes ambitious, when more 
than two electronic levels are taken into account. As ex- 
plained in Ref. HH, the number of paths that have to be 
considered is given by L 2 ( nc+1 >, where L is the number 
of electronic levels (i.e., here L = 4) and n c the number 
of retarded time steps that are accounted for. Obviously, 
keeping n c fixed, the number of paths is squared when 
one accounts for four instead of two electronic levels, and 
thus the numerical effort increases drastically. A reduc- 
tion of the number of paths can be achieved by invoking 
an on-the-fly path selection as described in Ref. HH and 
first introduced by Sim4^ This enables studies with a 
sufficiently fine time discretization but at the price of 
enforcing additional convergence tests, that have to be 
performed individually for different parameter ranges. 



C. Weak Coupling Theory 

In order to better understand the origin of the polar- 
ization dependence of the phonon-induced damping and 
some other features to be presented in Sec. IIII1 we set 
up equations of motion for the reduced electronic den- 
sity matrix within a weak coupling theory^. A more 
detailed description of this perturbative approach that is 
essentially equivalent to a standard Born-Markov master 
equation approach^ is given in the appendix. Here, we 
shall only give a short summary. 

First, we transform H as given in Eqs. (HJ-© from the 
electronic basis {|G), |+), |— ), |B)} to the basis spanned 
by the instantaneous eigenstates of H^ot+Hdot- light, that 
are usually referred to as dot-photon dressed states. To 
keep the formulas as simple as possible, we concentrate 
on the special case A = V cx = 0. For this choice, the 
dressed states are independent on f a± as well as on time 



and given by di = (1, 1, 1, 1), d 2 = (1, 1, -1, -1), d 3 = 
(1, —1, —1, 1) and d4 = (1, —1, 1, —1) with corresponding 
eigenenergies —h\j, where Ai = /+, A2 = /_, A3 = — /+, 
A 4 = -/_ and /± = ± U_) /2. 

Setting up equations of motion for the reduced elec- 
tronic density matrix pij = leads to a hierar- 
chy of phonon assisted quantities, that we truncate by 
factorizing at the level of two phonon assistances. Next, 
we integrate the one phonon assisted density matrix el- 
ements performing a Markov approximation and substi- 
tute the resulting expressions back into the equation for 
the reduced density matrix. Transforming back to the 
electronic basis {|G), |+), |— ), |B)} yields for example: 

d t pu = i/2 [/<t+ (P11-P22) + /<r_ (P14-P32)] (8) 
+ tt/4 [J(f a+ ) (pn+p 22 ) - Jtfa.) (Pu-Aa)] 
- tt/2 J(/ ct+ ) coth (HU + /2kT) Pl2 , 

where we have introduced the phonon spectral density 
q 

The main purpose for deriving the above approximate 
equations in addition to our numerically complete treat- 
ment within the path integral formalism is to obtain sim- 
ple closed form expressions for the relevant damping rates 
in terms of the phonon spectral density that allow for 
an explicit discussion of the dependencies on the field 
components f a± . The field strengths enter the damping 
rates because the energy conservation condition, that oc- 
curs in the Markov limit, enforces the phonon frequency 
to match the difference between electronic dressed state 
energies Ej = —K\j, which are directly related to f a± . 
Therefore, the phonon spectral density J(ui) is evalu- 
ated at the corresponding transition frequencies. For 
the special case A = V cx = discussed here, it is im- 
portant to note that J(w) is only evaluated at f a+ and 
f a _ [cf. Eq. (JSJ as an example], i.e., it is the strength of 
the cr_|_ and the er_ component separately, and not the 
total intensity, that matters. 

III. RESULTS 

In the following we will study the electronic dynamics 
that result from the interplay between optical excitation 
and the phonon-induced dephasing. For our calculations, 
that have been performed using the path integral ap- 
proach, we choose a e = 4nm, we set ah = 0.87a e , and 
assume that initially, the QD is in the unexcited ground 
state and the phonons in a thermal state at temperature 
T. A and V cx are treated as parameters. In Sec. MI AI 
we concentrate on linearly polarized excitations and in- 
vestigate how the laser pulse intensity affects the dynam- 
ics. Here, a special focus is put on comparing the results 
with those of the widely studied two-level model, where 
biexcitonic excitations are not accounted for. Sec. MI Bl 
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is devoted to the case of elliptically polarized excitations 
and the polarization dependence of the phonon influence. 



A. Linearly Polarized Excitations 

Let us first concentrate on linearly polarized excita- 
tions (f a+ = f a _ ) and consider the Rabi rotation scenario, 
where the final occupation after the pulse is recorded as 
a function of the applied pulse area a as given in Eq. (j7]l 
which is most convenient to be studied in experiments. 
Fig. [1] shows corresponding results for the ground state 
occupation Cg = Pii (green dashed line), the total single 
exciton occupation Ce = P22+P33 (red dotted line), and 
the biexciton occupation Cb = P44 (blue solid line) for a 
rectangular 12 ps lasting pulse at different temperatures 
for A = V cx = 0. In the phonon- free case [Fig. da)], the 
occupations of the ground and biexciton state perform 
undamped oscillations between and 1 with a period of 
2\/2ir, while the total single exciton occupation does not 
take values above 0.5 and oscillates twice as fast. Re- 
call, that the definition of the pulse area has been chosen 
such, that for a two level system the occupation oscil- 
lates with a period of 2tt and that for linearly polarized 
excitations f a± = //v2, explaining the additional factor 

\/2 that enters the period of the ground and biexciton 
state occupation. Accounting for the carrier-phonon in- 
teraction yields a damping and a renormalization of the 
Rabi period. As shown in Figs. mb)-[IJd) both effects 
strongly depend on temperature. Qualitatively similar 
results have previously been reported for electronic two- 
level systems22 1 24 1 28 1 46 

The probably most striking theoretical prediction for 
the dissipative two-level dynamics is that the phonon- 
induced damping should depend nonmonotonically on 
the pulse intensity resulting in a reappearance of Rabi 
rotations at high pulse area o 28 ' 38 . Recently, first experi- 
mental evidence for this reappearance has been reported^ 
for QDs driven by circularly polarized light, where a two- 
level model can be applied: the measured data showed 
a clear roll-off behaviour at high pulse areas. Of course, 
the question arises as to how this scenario appears for 
linear polarization, in particular for finite values of the 
biexciton binding energy, and which differences can be 
expected compared to the two-level case. To answer these 
questions, different Rabi rotation scenarios are presented 
in Fig.[5]for a 12 ps lasting rectangular pulse at T = 10 K. 
The rectangular pulse shape, that in experiments can be 
realized applying well known pulse-shaping techniques^, 
has been chosen as for such pulses the reappearance is 
more pronounced than for bell-shaped pulses^ 8 -. 

Shown in Fig.[^a) are Rabi rotations for A = V cyi = 0. 
While for small pulse areas the damping increases with 
rising pulse intensity (as also seen in Fig. [T]), it decreases 
at high pulse areas resulting in a reappearance of Rabi 
rotations. The pulse area a c with minimal Rabi rotation 
amplitude is the same for Cg, Ce and Cb and for the 
chosen parameters roughly given by a c = Hit. We stress 
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FIG. 1: (Color online) Occupation of the ground state Cg 
(green dashed line), the total single exciton occupation Ce 
(red dotted line), and the biexciton occupation Cb (blue solid 
line) as a function of the applied pulse area a for a linearly 
polarized rectangular pulse of 12 ps duration and A = 14 x = 0. 
(a) Displays results for the phonon-free case, (b)-(d) show the 
dynamics in the full model at temperatures of (b) T = 10, (c) 
30, and (d) 100 K. 



that the reappearance of Rabi rotations is different from 
the collapse and revival phenomenon within the Jaynes- 
Cummings model^, which describes the temporally pe- 
riodic carrier dynamics of a two-level system coupled to 
a single-mode quantized photon field. For the system 
considered here, the oscillation amplitudes of all occupa- 
tions decay monotonically with time, but with a damping 
depending nonmonotonically on the applied pulse area. 
The physical origin of this nonmonotonic damping, that 
is reflected in an undamping of Rabi rotations at high 
pulse areas, is the resonant nature of the carrier-phonon 
interaction^: the exciton-phonon coupling is maximal 
at intermediate pulse areas, where the frequencies of the 
electronic oscillations are resonant with the most strongly 
coupled lattice modes. 

In experimental studies, the final QD occupation after 
a pulse is often measured using a photocurrent detec- 
tion tcchniquc^ i 49 i 50 where under the action of an ap- 
plied electric field, the carriers, that have been created 
during the exciting pulse, tunnel from the QD result- 
ing in a photocurrent. For an exciton-biexciton system, 
this photocurrent is proportional to the weighted sum 
C := Ce + 2Cb of the generated single exciton and biex- 
citon populations. In the following, we will concentrate 
on this total electronic occupation instead of the single 
contributions Ce or Cb that in experiments are much 
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harder to extract individually. 

Plotted in Fig.^b) is C/2 for A = V" cx = (thick black 
line) together with the exciton occupation C*2ls of a two- 
level system (thin red line), that for V cx — we can 
model by choosing circularly polarized laser light. For 
the four-level case, the maximal signal is not only twice 
as high as in the two-level case (mind the scaling factor 
of 0.5 for C) but the pulse area with minimal ampli- 
tude, marking the beginning of the reappearance, is also 
considerably higher. This difference can be understood 
with the help of Eq. ([8|): the phonon- induced damping is 
expressed in terms of the phonon spectral density J{lo) 
that rises for small frequencies and peaks at a certain 
frequency uJ before it eventually decrease a 26 ' 38 . The im- 
portant point to note is that for the chosen parameters 
J is not evaluated at the total field strength / but at the 
circularly polarized field strengths f a+ and f a _ (which 
are identical for linear polarization). Further, for a given 
pulse area a, the circularly polarized field strengths of 
a linearly polarized pulse are by a factor of v2 smaller 
than the field strength of a circularly polarized pulse. 
Thus, to reach the field strength, where J and thereby 
the phonon-induced damping is maximal, higher pulse 
areas are needed for linearly polarized pulses than for 
circularly polarized pulses and hence, the reappearance 
starts later in the linear case. For the same reason, higher 
pulse areas are needed to fully restore the maximal Rabi 
rotation amplitude. 

Next, in Fig. [21c), we shall discuss the influence of the 
biexciton binding energy, which, depending on details of 
the quantum dot geometry, may take negative as well as 
positive values. While for A — (black dotted line) the 
dependence of C on a is similar to that of the two-level 
case, Rabi scenarios become much more complex for finite 
biexciton binding energies, where they qualitatively differ 
from the two- level results. For A = —1.0 (blue solid line) 
or +1.0 meV (blue dashed line) there is still a reappear- 
ance in the sense that first the Rabi rotation amplitude 
decreases with increasing pulse area before it eventually 
rises again, but in contrast to A = the signal consists 
of a superposition of different transition frequencies and 
the mean value varies with a. For small pulse areas, the 
mean value is less than one and the signal does almost 
not depend on the sign of A. For intermediate pulse ar- 
eas, it increases and we observe a pronounced difference 
between A = —1.0 and +1.0 meV, whereas for high pulse 
areas, both binding energies result again in similar oscil- 
lations with a mean value of approximately one. Further, 
there is a clearly visible increase of the oscillation period 
with rising pulse areas. 

To understand these features, it is helpful to consider 
the biexciton contribution 2Cb to the total photocurrent 
as shown in Fig. [5Jd). For small pulse areas, the driving 
field strengths f a± are much smaller than |A/7i|. From 
analytic expressions for driven two-level systems with- 
out a coupling to phonons 3 ^, it is well known that off- 
resonant transitions become unlikely, when the frequency 
detuning is larger than the driving field strength and that 
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FIG. 2: (Color online) Rabi rotation scenarios at T — 10 K for 
a linearly polarized rectangular pulse of 12 ps duration, (a) 
Cg (green dashed line), Ce (red dotted line), and Cb (blue 
solid line) for A = meV. (b) Total electronic occupation 
C — Ce + 2Cb (thick black line) and corresponding result for 
a two- level system (red thin line, see text), (c) C = Ce + 2Cb 
for A = (black dotted line), +1.0 (blue solid line), and 
— 1.0 meV (blue dashed line), (d) Biexciton contribution 2 Cb 
to the signals plotted in (c) for the same parameters as there. 



the dynamics does not depend on the sign of the detun- 
ing. Accounting for phonons. these statements remain 
approximatively valid as long as the time-scales consid- 
ered are much shorter than the time-scale, on which the 
system is driven to a stationary stated 3 -, which is fulfilled 
for small pulse areas. Thus, the biexciton is hardly ex- 
cited [cf. Fig. EJd)] in the regime of low pulse intensities 
and C basically given by Ce ■ With increasing pulse area, 
the generation of the biexciton becomes more and more 
likely and Cb starts to contribute, which is reflected in an 
increase of the mean value of C. At intermediate pulse 
areas, the contribution of Cb is noticeable but depends 
strongly on the sign of A. The latter can be explained as 
follows: as already mentioned above, the carrier-phonon 
interaction is maximal at intermediate pulse areas, where 
the system is quite rapidly driven to a stationary non- 
equilibrium state. As shown in Ref. l34l the stationary 
occupations depend on the question whether dark super- 
position states are realized or not and are also strongly 
dependent on A and T. For the given temperature of 
T = 10 K and a driving field strength of / = 2.0 ps -1 
[corresponding to pulse areas for which the Rabi rotation 
amplitude is minimal in Fig. 0c)], the stationary biexci- 
ton occupation is roughly 0.55 for A = 1 meV, while for 
A = —1.0 meV it is strikingly less and roughly given by 
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0.14. Although we never reach a stationary state within 
the considered pulse duration of 12 ps, fingerprints of 
the relaxation process become visible in the regime of 
an enhanced carrier-phonon interaction as given at inter- 
mediate pulse areas and lead to the strong dependence 
on the sign of A. At high pulse areas, the driving field 
strengths arc eventually significantly larger than |A//i|, 
the biexciton can be efficiently excited and due to the de- 
creasing phonon influence with increasing pulse area the 
system is more and more far from becoming stationary, 
which results in a diminishing dependence on the sign of 
A. Finally, the increasing period of the oscillations of 
C with rising a can be explained by the increasing con- 
tribution of Cb that shows a considerably longer period 
than Ce as already discussed in Fig. [T] and also shown in 
Fig. Ha). 



B. Elliptically Polarized Excitations 

So far, we have studied the exciton-biexciton dynam- 
ics under linearly polarized pulses. It was shown that 
the coupling between carriers and phonons strongly de- 
pends on the pulse intensity and that Rabi rotation sce- 
narios for finite A are qualitatively different from those 
of two-level models, while they are simliar for A ~ 0. 
In what follows, we shall demonstrate a remarkable po- 
larization dependence of the phonon-induced damping in 
an exciton-biexciton system, which has obviously no ana- 
logue in a two-level system. 

To quantify the ellipticity of the polarization, we first 
introduce a polarization parameter 7 via f a _ = jf a+ ■ For 
7 = 1, we obtain linearly polarized light, while a circu- 
larly polarized excitation is realized for 7 = 0. Plotted in 
Fig. Ela) is the total electronic occupation C = Ce + 2Cb 
as a function of time under cw-excitation for polariza- 
tions with 7 = 1 (black dashed line), 0.5 (red solid line), 
0.3 (green dotted line), and 0.1 (blue dottcd-dashed line), 
respectively. The total field strength has been fixed to 
/ = 2.0 ps -1 , i.e., the pulse intensity is the same for all 
displayed results, and we have chosen a vanishing biexci- 
ton binding energy A as well as a vanishing exchange in- 
teraction T4 X , as for A = V cx = an interpretation based 
on the simple rates derived in Sec. Ill Cl for this case can be 
made. For all excitation conditions, the total electronic 
occupation performs damped Rabi oscillations. But re- 
markably, the phonon-induced damping strongly depends 
on the polarization of the excitation. While the damping 
is strongest for linearly polarized light, the signal is the 
less damped the closer the polarization is to be circular. 
For 7 = 0.1, corresponding to an almost circularly po- 
larized excitation, the oscillation amplitude is after 50 ps 
more than 15 times larger than for linear polarization 
(7 = 1), where the oscillation amplitude is after 50 ps 
strongly reduced and the system already close to becom- 
ing stationary. As the carrier-phonon coupling as given 
by i?dot-ph is identical for er + and cr_ excitons, a polar- 
ization dependent damping of Rabi rotations is at first 
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FIG. 3: (Color online) Totally created charge C = C E + 2C B 
for different polarization parameters 7 (see text): 7 = 1.0 
(black dashed line), 0.5 (red solid line), 0.3 (green dotted 
line), and 0.1 (blue dashed-dotted line) at T = 10 K and for 
Vox = A = 0. (a) shows the dynamics under cw-excitation 
for a total field strength of / = 2.0 ps -1 as a function of time, 
(b) displays the final QD occupation after a Gaussian shaped 
pulse with a FWHM of 30 ps as a function of the applied pulse 
area. 



a very surprising finding. In light of a spin-independent 
exciton-phonon interaction one could rather expect that 
the damping is fully determined by the laser intensity and 
should thus depend only on / or other quantities related 
to the total intensity like f a+ + f a _ . However, our results 
clearly demonstrate that this is not the case and reveal 
a strong polarization dependence of the phonon-induced 
damping. 

For the discussion of this remarkable finding, cw- 
excitations are most advantageous, as for this choice the 
field strengths f a± are time-independent and an analy- 
sis based on field-strength dependent damping rates as 
given below is particularly simple. However, in experi- 
ments it is a challenge to obtain time-resolved traces of 
the electronic occupation as presented in Fig. Eta) and 
in most experimental studies pulsed excitations are used, 
where the final occupation is recorded as a function of 
the applied pulse area. We would like to stress, that the 
polarization dependence of the phonon-induced damping 
is also clearly seen for the latter scenario. This is ex- 
cmplarily shown in Fig.^b), where the final occupation 
is plotted as a function of the pulse area for a Gaussian 
shaped pulse with a pulse duration of 30 ps (FWHM): 
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the reduction of the Rabi rotation amplitude with rising 
pulse area is the less the closer the polarization is to the 
circular limit, reflecting a decreased damping in the time 
domain. 

The most convenient way to think of phonon-mediated 
relaxation is that the laser couples the excitonic states 
to form optically dressed states, that represent stable 
quasiparticles^. As discussed in Sec. Ill CI the damping 
of Rabi oscillations is determined by the phonon spectral 
density evaluated at the transition frequencies between 
the dressed states, which for A = V ex = are given by f a+ , 
/cr_j fa + — fa- an d /<t + +/<t_, but it turns out that terms 
involving f a+ + f tT _ or f a+ — f a _ drop out. This result 
explains the strong dependence of the damping on the 
polarization seen in Fig. [3] as follows (for simplicity, we 
assume again cw-excitation): while for linear polarization 
both circularly polarized components are given by // v2, 
for elliptical polarizations f a+ is larger and f a _ smaller 
than this value (without loss of generality we can assume 
f a+ to exceed corresponding to 7 < 1). Depending 
on the total laser intensity, we can then encounter two 
situations: either J{f a _ ) < J{fa + ) and the cr_ driven 
subsystem experiences a weaker damping than in the case 
of an equally strong linear driving [which is fulfilled for 
the driving strength of f — 2.0 ps -1 in Fig. |3Ja)], or 
J{fa+) < J(fa-) an d the damping of the <r + driven sub- 
system is less strong. In any case, one subsystem expe- 
riences a weaker damping than in the linear case, which 
results in a weaker damping of the total electronic occu- 
pation, to which both subsystems contribute. From this 
reasoning it follows straightforwardly that for clliptically 
polarized excitations the difference of the damping com- 
pared to the linear case should monotonically increase 
with rising ellipticity, which is indeed confirmed by the 
results shown in Fig. [3] 

For elliptical polarizations close to the circular limit, 
the carrier dynamics in one subsystem experiences even- 
tually hardly any damping and C oscillates with large 
amplitudes over a long time-interval. In Fig. EHa) this 
is nicely illustrated for 7 = 0.1 (blue dashed-dotted 
line). First, the signal rapidly oscillates. However, af- 
ter roughly 20 ps these fast oscillations have decayed. 
Subsequently, the dynamics changes qualitatively and is 
dominated by oscillations with a much longer period that 
are only weakly damped. While the fast oscillations re- 
flect the dynamics in the + driven subsystem, the much 
less damped oscillations with a longer period stem from 
the er_ driven subsystem. As for the chosen polarization 
fcr_ is roughly ten times smaller than f a+ , the oscillations 
in the er_ driven subsystem reveal a roughly ten times 
longer period and experience a considerably less damp- 
ing, because here J(/ CT _) <C J(/ CT+ ). This characteristic 
change in the dynamics is also expressed in corresponding 
results for pulsed excitations, cf. Fig. [3Jb). 

The limit of a strictly circular polarization with 7 = 
is special: for this excitation condition, one subsystem 
is driven with the total intensity, while the other sub- 
system is not driven at all. In consequence, while for 
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FIG. 4: (Color online) Totally created charge C = C E + 2C B 
under cw-excitation as a function of time at T — 10 K, / = 
2.0 ps -1 and V BX = for different polarizations with 7 = 1.0 
(black dashed line), 0.5 (red solid line), 0.3 (green dotted line), 
and 0.1 (blue dashed-dotted line) for situations, in which the 
relevant damping rates do not separatively depend on f a , or 
f a _ (see text). 



almost circularly polarized light the damping seen in the 
oscillations of the total occupation decreases with rising 
ellipticity and the time to become stationary rises with- 
out any limit for 7 — > 0, this time is finite for circularly 
polarized light as discussed in more detail in Figs. 2(b) 
and (c) of Ref. H 

It should be noted, that the finding that the damp- 
ing rate depends separately on f a+ and f„_ strictly 
holds only in the strong confinement limit as defined by 
Eqs. (HJ-© and in the limiting case A = V e * = dis- 
cussed so far. In the phonon-frec case, it is easily seen 
that for A = V cx = the model defined by Eqs. (P)-© can 
be decomposed into a product of two two- level systems, 
each evolving on its Bloch sphere with a Rabi frequency 
given by f a+ or f a _ , respectively. Thus, in this case we 
have two subsystems that evolve independently and re- 
act separately cither on the values of f a+ or f a _ . In a 
four-level system composed of two independent two-level 
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systems each coupled to phono ns, the phonon coupling 
of the two-exciton state must necessarily be twice the 
exciton-phonon coupling. Although this necessary con- 
dition is fulfilled in the exciton-biexciton system studied 
so far, we are actually not dealing with strictly indepen- 
dent dynamics on two Bloch-spheres as the two spheres 
arc coupled to the same phonon modes. The resulting in- 
terdependence of the dynamics on the two Bloch spheres 
can be seen, e.g., from the fact that for linearly polarized 
driving the stationary biexciton occupation reached at 
long times is not given by 1 /4 (which would be the prod- 
uct of the stationary values of 1/2 reached in any of the 
a± coupled two-level systems) but 1/3^ Nevertheless, 
even though a strict factorization into two independently 
evolving subsystems does not hold, a phonon-coupling 
of the biexciton state 7 q ,B of exactly twice the exciton 
phonon-coupling 7 q .± realizes a distinguished situation 
as in this case the phonon induced damping rates de- 
pend on f a+ and f a _ separately. However, this property 
is lost, when either 7 qj B ^ 27 q ,± or when the biexciton 
binding energy or the exchange interaction take finite 
values, or when the dipole moments associated with the 
ground state exciton transition Mq.± and the exciton 
biexciton transition M±.b are different. In each of these 
cases terms involving f a+ + f a _ or f a+ — f a _ enter the 
equations derived within the weak coupling theory. It 
turns out, however, that for many typical situations with 
deviations from the limiting case discussed so far, the 
amplitudes of the additional terms are small and the po- 
larization dependence is still close to the results shown in 
Fig. El This is illustrated in Fig. @] and shall be discussed 
in the remainder of the paper. 

In real systems, M±,b differs from Mg,± due to devia- 
tions from the strict strong confinement limit, where the 
multi-pair wave functions factorize in single particle wave 
functions for electrons and holes. Typical values observed 
in experiments indicate, that M±3 is roughly 5% smaller 
than Mq ; -|-i 31 i 51 For the same reason, 7 qj B ^ %Jq,± is 
possible. Figs. Hta) and BJb) show results for calcula- 
tions assuming M±,b = 0.95 Mg,± and 7 q ,B = 2.2 7 qj ±, 
respectively, for the same polarizations as in Fig. [3] It 
is clearly seen that the polarization dependence is only 
weakly influenced by these deviations from the strong 
confinement limit. 

Wc have also performed calculations (not shown) ac- 
counting for a finite exchange interaction of V cx = 
30 ficV, which represents a typical value. Although V ex 
yields a direct coupling of the two subspaces driven by a± 
polarized light, we found only very small changes com- 
pared to V ex = 0, except for strict circular polarization^, 
or in the regime of nearly circularly polarized excitations, 
when the time-scale on which the total system is driven 
to a stationary state increases drastically and exceeds 
h/V„. 

Finally, let us turn to the influence of finite values of 
the biexciton binding energy, as illustrated in Figs. @Jc) 
and HJd) f° r A = 1.0 and —1.0 mcV, respectively. At 
first one may notice that the mean value of the oscilla- 



tions decreases with rising ellipticity of the polarization. 
This lowering of the total signal is due to the fact, that 
the biexciton generation becomes more unlikely, when 
/cr_ decreases, as in this case the exciton biexciton tran- 
sition given an exciton generated by er + polarized light 
becomes less probable. In addition, there is no longer a 
cross-over between fast and slow oscillations for almost 
circularly polarized excitations. However, despite these 
differences compared to the case of a vanishing biexciton 
binding energy, we find a very similar behavior as dis- 
cussed so far: the decay of the oscillation amplitude is 
the less pronounced the more far the excitation is from 
the linear case, clearly illustrating a strong dependence of 
the phonon-induccd damping on the polarization of the 
driving laser light. This picture changes, however, for 
larger A, which are still feasible, as shown in Fig. @Jc) 
for A = 5 mcV. Here, the polarization dependence is 
drastically suppressed. As the coupling to the biexciton 
decreases with rising A, this indicates, that an efficient 
coupling to the biexciton is of fundamental importance 
for the polarization dependence of the damping. Indeed, 
when we completely decouple the biexciton from the dy- 
namics by setting by hand M±.b = 0, then the polar- 
ization dependence vanishes altogether and the phonon- 
induced damping is fully determined by the pulse inten- 
sity [cf. Fig. |U(f), where all four lines coincide]. Only 
when the biexciton is coupled, there is the possibility to 
essentially drive two different Bloch spheres with individ- 
ual Rabi frequencies within the four-level system consid- 
ered here, giving rise to the polarization dependence of 
the phonon-induced damping as explained above. 



IV. CONCLUSIONS 

In summary, we have presented a numerically com- 
plete analysis of the phonon impact on optically driven 
exciton-biexciton QD systems. Rabi rotation scenarios 
for finite biexciton binding energies were shown to dif- 
fer significantly from the two-level case and offer mani- 
fold insights in the dynamics of an exciton-biexciton sys- 
tem, including the interplay between off-resonant driv- 
ing and quantum dissipative relaxation. In addition, we 
have demonstrated that although the carrier-phonon in- 
teraction is independent of the spin, the phonon-induced 
damping does not only depend on the pulse intensity but 
is for typical parameters also strongly affected by the po- 
larization of the driving laser pulse. While the damping 
is strongest for linearly polarized excitations, it is the 
more reduced the closer the polarization is to the cir- 
cular limit. For almost circularly polarized excitations 
and small biexciton binding energies, the Rabi scenario 
changes strongly in the course of time reflecting differ- 
ent dynamics in the a± driven subsystems. The origin 
of the polarization dependent phonon-induced damping 
is the possibility to drive two Bloch spheres essentially 
each with its own Rabi-frequcncy f a± that arc deter- 
mined by the polarization, together with the fact that 
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it is the Rabi frequency that determines the damping 
in a two-level system. This possibility is intimately re- 
lated to an efficient coupling to the biexciton state, which 
highlights the pivotal role of this two-pair state for the 
discussed dynamics. 
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VI. APPENDIX 

In this appendix, we give a short description of the 
weak coupling theory presented in Sec. Ill CI concentrat- 
ing on the limit of strong electronic confinement as de- 
fined by Eqs. CO)-© and the special case A = y C x = 0. 
The Heisenberg equations of motion for the electronic 
operators p,^ = \di)(dj \ read 

dtPij = (10) 

where the overbar indicates, that we are working in the 
dressed state basis. By taking expectation values we ob- 
tain equations of motion for the elements of the reduced 
density matrix Ji^ — (p^ ) , where single phonon assisted 
density matrix elements of the form {Pijb^) or (Pijbq) 
enter. Setting up equations of motion for these single 
phonon assisted quantities, in turn, introduces double 
phonon assistances. In order to truncate the resulting 
hierarchy of higher order phonon assisted quantities, we 
set 

(fijtfqbq') = Pij (b^bq') = pij <5 q , q ' (11a) 
(^A 6 q'> = Pij (Mq') = Pij K + 1) <W (Hb) 
&AV> = (^4 & q'>= ' ( UC ) 

where n q = l/(exp(fiw q /fcBT) — 1) denotes the Bosc- 
Einstcin occupancy at temperature T. By doing so, one 
obtains for the single phonon assisted elements a set of 
closed equations. Exemplarily, we have 

9t(Pi 2 bq) =i(-w q + A 2 i)(p" 12 V (12) 

+ -Ifi-VPn - "q(P22 +P42) 
+ ("q + 1 )(Pll +Pl3)] ■ 

Integrating these equations yields integrals of the form 

/>oo 

/ exp(iwT)p i:j (t - T)dr . (13) 
Jo 

Next, we perform a Markov approximation by identifying 
slowly varying variables and taking them out of these 



memory integrals. While occupations can be regarded as 
slowly varying, this is usually not the case for coherences, 
that are generally rapidly changing quantities. However, 
as their free oscillation frequency is given by Xji := Xj — 
Ai, we can set ~p~ij(t) = cx ] p(i\jit)'p i j{t), where ~p~ij(t) can 
be expected to vary slowly. Taking these variables out of 
the integral yields 

/>oo />oo 

/ exp(iwT)pj J (t— r)dr = ~p i3 (t) / exp[i(w — Xji)r]dT 
Jo Jo 

(14) 

and the remaining integrals can be easily evaluated using 

r 00 p 

/ exp(iu)T)dT = it5(lj) + i — , (15) 
Jo ^ 

where P denotes Cauchy's principal value. It is well 
known that the principal value terms only lead to a small 
rcnormalization of the system Hamiltonian but do not 
directly contribute to dephasing.— Therefore, we neglect 
these terms and obtain: 

(Pl2&q) = -^r\-~ 2 Pl2^(^q) - "-qP22^( W q + ^1 2 ) (16) 

- ?l qP42^( w q + ^u) 

+ (71 q + l)p n 5(w q + Al 2 ) 

+ («q + l)Pl3^( W q ^ A 23)] ■ 

Here, only A i4 = A 23 = U+ and \ 12 = A 43 = f a _ ap- 
pear, although in Eq. (fT2j) also combinations like A24 = 
/<j + — f a _ (via the free oscillation of (P42)) entered, be- 
cause the corresponding contributions cancel out. It is 
an important point to note, that this holds for all other 
phonon-assisted density matrix elements as well: from 
here on only 8{uj — f a+ ) and 5(lo — f rT _ ) enter the re- 
sulting equations. Terms involving 8[uj — {f a+ + /<r_)] 
or 5[uj — (f a+ — f a _)\ are absent due to the structure of 
#dot-iight and i?dot-phonon in the special case considered 
here. Introducing the phonon spectral density 

■7M -Wq) ( 17 ) 

q 

and substituting the expressions for the single phonon 
assisted density matrix elements into Eq. (|10p eventu- 
ally leads to a closed set of equations of motion for the 
elements of the reduced electronic density matrix. 

We mention, that the weak coupling theory derived 
here yields the same result as a four-level extension of 
the master-equation approach that has been presented in 
Ref. H for the case of an electronic two-level system. For 
two-level systems the latter approach has been compared 
to a much more advanced variational master-equation 
approach in Ref. 0, which in turn has been shown to 
agree well with numerically exact path integral calcula- 
tions for not too high temperatures and standard GaAs 
parameters, that represent a weak carrier-phonon cou- 
pling regime. We checked that in the four-level case, the 
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deviations of the weak coupling theory to numerically 
exact calculations are qualitatively similar as in the two- 
level case and refer to Ref. @ for a comparison. Readers 
interested in the validity ranges of approximate meth- 
ods in the regimes of strong carrier-phonon couplings 
or high temperatures are referred to Refs. and HH, 



where path integral calculations have been compared to 
a second-order and fourth-order correlation expansion, 
clearly demonstrating the superiority of a numerically 
exact approach, accounting fully for all non-Markovian 
effects and arbitrary multiphonon processes. 
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